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Abstract. We present a simultaneous Chandra and XMM-Newton observation of the Castor sextett, focusing on Castor A and 
Castor B, two spectroscopic binaries with early-type primaries. Our study represents the first unambiguous X-ray detection of 
all three visual components in the Castor sextett making use of the unprecedented spatial resolution of Chandra. Of the present 
day X-ray instruments only Chandra can isolate the X-ray lightcurves and spectra of Castor A and B (angular separation ~ 4"). 
We compare the Chandra observation with XMM-Newton data obtained simultaneously. Albeit not able to resolve Castor A and 
Castor B from each other, the higher sensitivity of XMM-Newton allows for a quantitative analysis of their combined high- 
resolution spectrum. He-like line triplets are used to examine the temperature and the density in the corona of Castor AB. The 
oxygen triplet provides a density of ric — (0.5...1) 10^^° cm~^, typical for stellar coronae. The analysis for the neon triplet 
results in much higher densities. By means of a simulated RGS spectrum we estimate the contaminating effect of iron lines to 
the neon triplet. The temporal variability of Castor AB is studied using data collected with the European Photon Imaging Camera 
onboard XMM-Newton. Strong flare activity is observed with typical rise times of ~ 10 min and exponential decays which are 
by a factor of 2 — 3 slower than the rise. Combining the data acquired simultaneously with Chandra and XMM-Newton each flare 
can be assigned to its host. Thus we verify that both Castor A and Castor B exhibit flares. Our comparison with the conditions 
of the coronal plasma of other stars shows that Castor AB behave like typical late-type coronal X-ray emitters supporting the 
common notion that the late-type secondaries within each spectroscopic binary are the sites of the X-ray production. 
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^ 1. Introduction have argued in favor of Castor A as the X-ray source be- 

• • cause they detected microwaves from Castor A but not from 

_ > Castor is a sextett system at a distance of 16 pc, and composed Castor B. The latter was, however, subsequently detected as a 

^ , of three visual stars each of which is by itself a spectroscopic ^^^.^ ^^^^^^ j Oiidel & Schmitt 199 6l. X-ray and mi- 

^ . binary. Castor C (= YY Gem) consists of two late-type stars ^^^^^^^ ^^^^.^^ ^^p^^^^^ associated with each other 

Ci ■ (dMle + dMle) and is separated from Castor AB by ^ 74 . ^^^^^^^ ^^^j^ generated as a result of magnetic energy re- 

The primaries in both Castor A and Castor B are of spectral j^^^^ ^^^.^^ ^^^^^^^^ ^^^^ ^^^^^ ^^^^^^ ^_ 

type A (Al V and A5 V, respectively). Neither strong stellar ^^.^^^^^ ^^^^^^.^ ^1^^ ^^^-^^ pj^^^^^ quiescence. X-ray and 

winds (typical for early-type stars) nor magnetic activity (typ- i^icrowave emission in stars have been found to be strongly 

ical for late-type stars) are thought to produce X-ray emission ^^^related jBenz & Giidel 19 94..Gudel & Benz 1993 1 . 
in early A-type stars. Nevertheless, observations by EXOSAT 

allowed to separate two sources of X-ray emission, YYGem I" ^ recent XMM-Newton observation [Hlidel et al. (200^1 

and the visual binary of Castor AB (Pallavicini et al. 1990ai. ^und *at all three visual binaries in the sextett are X-ray 

Castor AB was also margin ally resolved from YYGem by sources. The angular separation of Castor A andB is - 4", 

ASCA (see lGotthelt'et^rW94i . The most likely explanation ^t the limit of the spatial resolving power of XMM-Newton. 

is that the presumably late-type spectroscopic companions The superior resolution of Chandra now allows us to present 

of Castor A and/or Castor B are responsible for the X-ray a" unambiguous X-ray detection of both Castor A and 

emission. The companion of Castor A is most likely a late- Castor B, and for the first time to isolate the X-ray spec- 

K star in a 9.21 d eccentric orbit, while CastorB's compan- ^rum and lightcurves of both binaries. A contemporaneous 

ion seems to be an early-M dwarf with a 2.93 d circular or- XMM-Newton observation provides higher signal-to-noise 

bit dOiidel & ^chmitt l996limintz l^}j}jt . |53Fmitt et al. (1^94)| at both medium- and high-spectral resolution. To date high- 

resolution X-ray spectra have been presented for about a 

Send offprint requests to: B. Stelzer dozen late-type Stars (e.g. Giidel et al. 2001 Ness et al. 2001' 

Correspondence to: B. Stelzer, stelzer@astropa.unipa.it IMewe et al. 20011 iRaassen et al. 2002. .Stelzer et al. 2002. 
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INess et al. 2002at . providing for the first time direct access to 
the physical conditions in stellar coronae by measuring fluxes 
of individual emission lines. The relative strength of lines from 
different elements shed light on the temperature, density and 
abundances in the emitting region. Enlarging the data base 
of high-resolution studies of stars is important to distinguish 
which of the observed properties are typical for stellar coronae 
and which of the objects are peculiar Observing a variety of 
stars with different spectral type, multiplicity, rotation, and 
age will eventually allow to disentangle the influences of 
these parameters on stellar activity. The combination of our 
simultaneous XMM-Newton and Chandra observation provides 
a unique data set with all the capabilities of the present-day 
generation of X-ray instruments, and enables a detailed study 
of the characteristics of three coronal X-ray sources in the 
Castor system. 

In Sect. |2] we introduce the observations. The data reduc- 
tion is described in Sect. |3] We present in brief our analysis 
of the medium-resolution spectrum (Sect. 1^. Modelling the 
high-resolution spectra obtained with both Chandra and XMM- 
Newton, and interpretation of the results is the major purpose 
of this paper (Sect. |5|l- Therefore, Sect. |4] does not aim at an 
exhaustive discussion of the medium-resolution spectrum of 
Castor AB. Rather it is introduced because of the additional 
information it provides for a better understanding of the high- 
resolution spectrum which is used along this paper A study 
of the X-ray variability of Castor A and B based on the X-ray 
lightcurves is given in Sect. |6l In Sect. we compare the ac- 
tivity of Castor A and B to other active stars, and in Sect. Owe 
summarize our results. 

2. Observations 

We have obtained deep X-ray observations of the Castor sys- 
tem with both Chandra and XMM-Newton (Obs-IDs 28 and 
0112880801, respectively). The observations took place si- 
multaneously on Sep 29/30, 2000. Chandra was used in the 
the LETGS configuration, i.e. the Low Energy Transmission 
Grating (LETG) was combined with the High Resolution 
Camera for Spectroscopy (HRC-S). XMM-Newton allows 
to perform high-resolution spectroscopy with the Reflection 
Grating Spectrometer (RGS) and at the same time CCD imag- 
ing and spectroscopy with the European Photon Imaging 
Camera (EPIC). During the observation EPIC was operated in 
full-window mode with the thick filter inserted due to the op- 
tical brightness of Castor {V ^ 1.6 mag). For the same reason 
the Optical Monitor was in closed position. 

The prime target of this campaign was YYGem (= 
Castor C), discussed by |Stelzer et al. (2002)| where we also pre- 
sented the observing log. The instrumental setup allows to anal- 
yse also the temporal and spectral behavior of the other two vi- 
sual components of the sextett. Castor A and Castor B. In this 
paper we focus on the X-ray emission of these two binary stars. 

3. Data Reduction 

Fig-fflshows the zeroth order LETGS image of Castor A and B. 
Their separation on the LETGS image is 3.8" and the po- 
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Fig. 1. Chandra zeroth order LETGS image of Castor A and 
Castor B. Their separation in this image is 3.8". 



sition angle is 64.9°, in agreement with the value expected 
from the optical ephemeris lIHeintz 1988t . With the LETGS 
both the dispersed spectrum of Castor A and Castor B and their 
lightcurves are separable. In Fig. |2] the raw LETGS spectrum 
of the Castor system is displayed. To extract the LETGS spec- 
trum for Castor A and Castor B we followed the CIAO 2.2 sci- 
ence threads' for grating spectroscopy with multiple sources. 
Three sources are detected with tgdetect, corresponding to 
YYGem, Castor A, and Castor B. We selected a rectangular 
box of 0.0005° width in the cross-dispersion direction defin- 
ing the source spectrum of each Castor A and Castor B. The 
background was extracted from two larger rectangular regions 
above and below the source spectrum, and excluding respec- 
tively the other two X-ray sources. Lightcurves for Castor A 
and Castor B are extracted from a circle of radius 2" at their ze- 
roth order position. The background is estimated from an aver- 
age of the counts in four square-shaped regions arranged sym- 
metrically around the source. After background subtraction we 
find a total of 1989±46 cts and 1667±43 cts for Castor A and B 
respectively. The exposure time of Chandra was 57.9 ksec. 

With XMM-Newton Castor A and B are not resolved, but 
its instruments provide higher sensitivity allowing for a de- 
tailed and quantitative analysis. The RGS data were reduced 
with the standard pipeline rgsproc of the XMM-Newton Science 
Analysis System (SAS), version 5.3. With help of the bkgex- 
clude parameter we made sure that YY Gem is excluded from 
the background extraction region. We extracted the first order 
total spectrum of Castor AB (source plus background) and a 
background spectrum for analysis with the CORA ^ line fitting 
program (see Sect.|5jl. The RGS response matrix was generated 
with 6500 energy bins. 

' http://asc.harvard.edu/ciao/threads/spectra-multi-hrc/ 
" CORA is provided by J.-U. Ness and available at 
http://www.hs.uni-hamburg.de/DE/Ins/Per/Ness/Cora/index.html_ 
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Fig. 2. Raw LETGS spectrum of the Castor system. 



The data collected by the European Photon Imaging 
Camera (EPIC) onboard XMM-Newton was also analysed with 
the standard ZMM-SAS tools, version 5.3. The position of 
the X-ray source representing Castor AB was determined by 
searching in an iterative way for the photon centroid at the ex- 
pected position. Lightcurves and spectra for Castor AB were 
extracted from a 25" sized circle centered on their X-ray po- 
sition. This value for the extraction radius corresponds to the 
point where the integrated radial distribution of counts around 
the position of Castor AB flattens out. This way we ensure that 
YY Gem does not contribute significantly to the selected pho- 
tons. A ~ 1 ksec long time interval near the end of the XMM- 
Newton pointing was cut out because of higher background (see 
Fig.|6|l. The total useful exposure time with EPIC pn amounts 
to 48.6 ksec. 



Table 1. Best fit parameters of a 3-temperature VMEKAL 
model describing the EPIC pn spectrum of Castor AB. 



xlAdof) m m m igEUi 

[keV] [keV] [keV] [cm^^] 

1.15(525) 0.27^ni 0.78to°ot '^■^^-oll ^^-^^-ox. 

IgEM'z IgEMs O Ne Fe 

OO+0.05 nr,+0.05 n cq+O OS ^^.4+0.22 () M+O ^ 



4. Medium-Resolution Spectrum: EPIC pn 

Since this paper concentrates on the specific new achievements 
enabled by the high spatial and spectral resolution of Chandra 
and XMM-Newton the EPIC spectrum will only briefly be in- 
troduced. We discuss only the pn spectrum, which provides the 
highest sensitivity among the EPIC instruments. 

For the spectral fitting we select photons from a circu- 
lar region as described in Sect. |3] A background spectrum is 
obtained from a nearby source-free region of the same area. 
We restrict the analysis to energies above O.SkeV, because the 
spectral response of pn is not well understood for lower ener- 
gies. We use the pn redistribution matrix for single and double 
event patterns released in April 2002 together with an ancil- 
liary response file generated in the course of our data reduction 
process. Modelling of the spectrum is performed with XSPEC 
version 11. 2.0. 



Considerable emission is seen in the pn spectrum up to 
comparatively high energies. The statistics are good enough 
to visually identify an Fe feature at 6.7 keV. The shape of 
the spectrum is best approximated by a three-temperature (3- 
T) model describing thermal emission from a hot, optically thin 
plasma (three VMEKAL models). In the selected energy range 
absorption is negligible. A model with solar abundances gives 
Xrcd = ^-^^ f^i" 528 dof, and leaves substantial residuals espe- 
cially near 10 A(neon) and 22A(oxygen). An acceptable solu- 
tion (Xrod ^ ^) found for subsolar abundances of iron, oxy- 
gen, and neon. All other elements do not have strong enough 
lines in the sensitive spectral range of EPIC to change the fit 
result significantly. We summarize the best fit parameters for 
the 3-T model with variable abundances in Tabled 
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5. High-Resolution Spectra: LETGS and RGS 

The high-resolution spectrum of both Castor A and Castor B 
(displayed in Fig.|3} is that of a typical active late-type dwarf 
star The LETGS spectrum, though being quite weak, shows 
that the major emission lines in both stars are from OVIII, 
OVII, NelX, and FeXVII. A quantitative analysis of the 
LETGS spectrum by measuring fluxes for individual lines is 
subject to large uncertainties because the effective area of the 
LETGS for off-axis sources is not well-understood. We esti- 
mate these uncertainties (a sum of HRC-S and grating effi- 
ciency, transmission in the spatially inhomogeneous filter) to 
amount to ^ 10 %. This is comparable to or smaller than the 
statistical errors due to the faintness of both Castor binaries. 
RGS provides higher sensitivity, but no information on the lo- 
cation of the X-ray emission within the Castor AB system. In 
view of the similarity of the LETGS spectrum of Castor A and 
Castor B, we feel justified to model the combined RGS spec- 
trum, and consider it as representative for both stars. 

In the following we concentrate on the lines from 
hydrogen- and helium-like ions. The triplets of helium-like ions 
are composed of resonance r, intercombination i, and forbid- 
den / Une, and represent the most interesting lines in terms 
of plasma diagnostics: Their line ratios G — {f + i)/r and 
R — //i are sensitive to the electron temperature To and 
the electron density of the emitting plasma, respectively 
JGabriel & Jordan 1969> . The ratio of hydrogen-like Lya and 
helium-like r line provides a measure for the coronal tempera- 
ture, as it reflects the presence of different ionization stages of 
the same element. In the case of Castor AB the only two ele- 
ments with considerable emission above the background in the 
triplets and Lya lines are oxygen and neon. The data analysis 
and results for these two elements are discussed below. 

We measured the line strengths with the CORA line fit- 
ting application (Ness et al. 2001 1. For the LETGS data we 
used Gaussians to represent the individual emission lines. For 
RGS line shapes Lorentzians are a better representation than 
Gaussians. All lines are added to a constant describing the lo- 
cal background. 

5. 1. Oxygen 

In Fig. 13 we show the LETGS and the RGS spectrum in the 
region of the O VII triplet. 

To increase the signal both sides of the LETGS spectrum 
were co-added. To guide the fit we held the distance between 
the Une centers fixed on the expected value, and restricted the 
line widths to 0.03 A corresponding to the instrumental width 
of the LETGS. This latter restriction assures that random sta- 
tistical fluctuations do not influence the line profile. Given the 
large uncertainties related to poor statistics and low effective 
area at the off-axis position a conversion of the LETGS counts 
to line fluxes would not warrant meaningful results, and we re- 
nounce on performing this exercise. 

For the quantitative analysis only the RGS is used. As the 
CCD in the relevant part of the spectrum has failed on RGS 2, 
only RGS 1 is available for the analysis of the O VII triplet. Due 
to the higher quality of the RGS spectrum as compared to the 





Fig. 4. Triplet of helium-like oxygen. From top to bottom: 
Castor A with LETGS, Castor B with LETGS, Castor AB with 
RGS. Note the different bin size. 
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Fig. 3. First order RGS count rate spectrum (top) of Castor AB and LETGS count rate spectra of Castor A and B (middle and 
bottom). The counts from RGS 1 and RGS 2 have been added. Note, that due to the failed CCDs on each RGS the relative 
strengths of some lines appear distorted in the combined spectrum. For the LETGS left and right side have been added. Outside 
the displayed range only one emission line is detected, the Lya line of C VI at 33.7 A. 



LETGS it was possible to leave the line centers and line widths 
free in the fitting process. In Table|2we give the line parameters 
(central wavelength, width, number of counts, and photon flux) 
for the OVII triplet and OVIII Lya line. The flux has been 
determined from the effective area of the RGS at the respective 
wavelength after modeling of the count spectrum with CORA. 



5.2. Neon 

The neon triplet is more difficult to model, as it is known to be 
blended with a number of iron lines. We included additional 
Lorentzians to account for iron emission at 13.79 A (from 
Fe XIX) and 13.83 A (from Fe XVII). Further lines of Fe XIX 
may contaminate the r and the i line of Ne IX, but are indistin- 
guishable from the latter ones in the data. As a first approxima- 
tion we performed the formal analysis of the neon triplet in the 
same way as for oxygen, assuming that the contamination by 
iron is negligible. The best fit together with the data is shown in 
Fig-Elon the top. Again only one RGS instrument is available, 
RGS 2, because of a CCD failure in the relevant spectral range 
on RGS 1. 



To investigate the effect of iron contamination in the neon 
triplet we made use of the information obtained from the 
medium-resolution EPIC spectrum (discussed in more detail 
in Sect. |4}. Modeling the EPIC spectrum provides an esti- 
mate for temperature and abundances in the X-ray emitting re- 
gion, but the density can not be well constrained. The contri- 
bution of iron to the line emission in the critical spectral range 
(A - 13.4.. 13. 8A) can be estimated by generating an artifi- 
cial RGS spectrum for a plasma that has the same temperature 
structure as the EPIC spectrum, but no neon. We performed 
this simulation with XSPEC setting the Ne abundance of the 
EPIC 3-T best fit model of Table^to zero. Then we subtracted 
the simulated neon-free spectrum from the data. The result is 
shown in Fig. |5] (bottom), and should give an idea about the 
actual strength of the Ne IX triplet. Clearly the r, i, and / lines 
are more pronounced as compared to the pure data, implying 
that the iron contamination of the data is substantial render- 
ing any physical parameters derived from the spectrum at these 
wavelengths questionable. 
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Table 2. Line parameters for the Lya line of the H-like ion and 
triplet lines of the He-like ion of oxygen. The best fit result 
from both RGS and LETGS data are listed. Fluxes were com- 
puted only for RGS, because of the uncertain effective area of 
the LETGS and low signal-to-noise in the spectrum. 



Identif. 



Line Parameters 





A 


(7 


/ 


Ph.flux [10"^ 




[A] 


[A] 


Fctsl 


cts /s/ cm^ ] 




Castor AB 


- RGS 1 




O VTTT T vrv 

V.y V ill -1— 'V Lt 


18.973 


0.087 


71 7 ^ + 9Q 4 


oc 9 _|_ 1 n 


OVIIr 


21.605 


0.085 


218.0 ± 17.7 


8.8 ± 0.7 


OVIIi 


21.811 


0.048 


50.0 ± 10.0 


2.1 ±0.4 


OVII/ 


22.091 


0.071 


138.0 ± 14.3 


5.7 ±0.6 




Castor A - 


LETGS 




VIII Lya 


18.971 


0.033 


170.2 ± 14.1 




OVIIr 


21.608 


0.030 


44.2 ± 8.1 




OVIH 


= 21.808 


0.027 


15.5 ±5.8 




OVII/ 


= 22.108 


0.027 


30.1 ± 7.2 






Castor B - 


LETGS 




O VIII Lya 


18.963 


0.025 


127.3 ± 12.2 




OVIIr 


21.587 


0.031 


27.8 ±7.1 




OVIIi 


= 21.787 


0.030 


17.3 ± 6.2 




OVII/ 


= 22.087 


0.027 


26.7 ±6.7 





5.3. Coronal Density 

We derive the electron temperature To and density ric for 
the oxygen and neon triplets using the calculations by 
[Porquet et al. (20(JT}] for a collision dominated plasma. As dis- 
cussed above the Une fluxes for the neon triplet derived from 
the data are severely affected by contaminating iron lines. We 
continue to include Ne IX in the analysis merely with the aim 
to investigate the effect of line blending in this spectral range. 

The assumption of a coUisional plasma is justified if G ~ 1, 
while larger values of G prevail in photo-ionized plasmas. UV 
radiation from the star can modify the relative intensities of 
the emission lines by populating one level at the expense of 
another one (see Porquet et al. 2001,) . In particular, the upper 
level of the forbidden transition can be depopulated to the up- 
per level of the intercombination transition. For typical late- 
type coronal X-ray emitters the UV flux is low and can be ne- 
glected. However, for Castor A and Castor B the A-type pri- 
maries in each of the binaries provide substantial UV flux at 
the position of the X-ray emitting secondary. The influence of 
radiative processes can be expressed in the radiation tempera- 
ture Trad- We compute Tiad making use of UV fluxes derived 
from archived International Ultraviolet Explorer (lUE) data of 
Castor A and B. In Table|3]we tabulate the lUE fluxes of each of 
the two binaries at the wavelengths corresponding to the excita- 
tion energy between the upper levels of the / and the i line for 
oxygen and neon. Making use of the distance (16 pc), stellar ra- 
dius (derived from the Stefan-Boltzmann law), and limb dark- 
ening (from Table 1 of .Diaz-Cordoves et al. 19951 of Castor A 
and B we convert the UV flux to intensity and derive Trad from 
the Planck curve for a blackbody emitting at the corresponding 
UV wavelength. The effect of the UV emission at the position 




Fig. 5. Triplet of helium-like neon, top - Castor AB observed 
with RGS, fit includes two iron lines that are blended with 
the triplet (see text in Sect. 15. 2> : bottom - RGS spectrum of 
Castor AB after subtraction of a simulated neon-free spectrum, 
i.e. after removal of the contribution of iron. 



of the secondary is proportionate to the 'dilution factor' W , 
which depends on the stellar radius and the distance to the UV 
source, i.e. the binary separation (see [Porquet et al. 2001] for 
the definition of W). 

The case of Castor is complicated by the fact that the UV 
field and geometrical situation of both components (A and 
B) differ slightly from each other, but only a combined mea- 
surement of the X-ray line ratios is available. Our RGS mea- 
surements for the G, R, and Lya/r ratios are summarized 
in Table |3 For the G ratio the effect of UV photo-excitation 
is negligible and we can infer an average temperature for 
Castor AB. We find ~ 2 MK for oxygen and ~ 8 MK for neon. 
Using these temperatures we calculated the R ratio as a func- 
tion of density for the two components separately: (A) adopt- 
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Table 3. Stellar parameters, lUE flux and radiation tempera- 
ture of Castor A and B at the wavelengths corresponding to the 
/ — > i line transition which is critical for the determination 
of coronal densities (see text in Sect. 15. 3> . 



Castor A 



Castor B 



effective temperature (1) 

gravity (1) 

radius 

limb darkening coeff. (2) 
binary separation (3) 
dilution factor 



Teff = 10286 K 

lgg = 4.2 
i?. = 2.4 i?0 
e = 0.5 
aAa = 0.121 AU 
W = 0.002 



Teff = 8842 K 

lg3 = 4.0 
Ft* =3.3 7?0 
e = 0.5 
aBb = 0.052 AU 
W = 0.02 



sity have not changed by much. However, the flux of i and /, 
and subsequently any physical parameters derived from them, 
are subject to large uncertainties. To summarize, it can be con- 
cluded that physical parameters derived from the Ne IX triplet 
should always be regarded with caution, and this part of the 
spectrum is not suited for an investigation of the conditions in 
the corona with the RGS instrument. Among the current X-ray 
instruments only the HETG is able to sufficiently resolve a stel- 
lar spectrum in this region (see lNess et al. 2002b> . 

6. X-ray Lightcurves 



OVII 



NelX OVII 



NelX 



Afl [A] 1637 1270 1637 1270 

/uv [10^°erg/s/cmVA] 8 4.5 8.5 4 

r,ad [K] 9950 8730 9280 8250 

(1) -ICayrel de Str obel et al. 1980| (2) -|Diaz-Cord oves et al. 19951 
(3) - M. Giidel, priv. comm., * V. Costa, priv. comm. 



Table 4. Line ratios derived from the unresolved RGS spec- 
trum of Castor AB and the RGS spectrum cleaned from iron- 
contribution (see text): G = {f + i)/r, R = f/i, and 
Lya/r. Corresponding electron temperature To and electron 
density ric are derived using the calculations for CIE by 
Iforquetet al. {imT]\ 



G 


R 


Lya/r 


T 


Tie 








[MK] 


[lO^^cm-^] 


Oxygen 


0.88 ±0.11 


2.79 ±0.63 


2.86 ±0.26 


~ 2 


(0.5... 1) 


Neon 


0.43 ± 0.09 


1.55 ±0.61 


1.04 ±0.13 


~ 8 


100 




Neon, 


data - simulation 






0.54 ±0.16 


1.10 ±0.56 




~ 7 


200 



ing the radiation temperature and dilution factor of Castor A, 
and (B) for the values of T^ad and W of Castor B. While for 
neon Tj-ad < 10000 K has no effect on the R ratio, for oxygen 
there are slight differences between the two calculations: We 



find rip 



10 cm for Castor A, and = 5 10 cm for 



Castor B. Note that this result was derived under the assump- 
tion that both stars exhibit the same line ratios. Slight differ- 
ences in the line ratios may lead to an effect on the density of 
the same order as the range given above. We summarize the re- 
sults for the temperatures and densities in Table0]together with 
the line ratios. 

The RGS data of the Ne IX triplet results in a very high 
density. Such high densities are untypical, though not gener- 
ally excluded for active stars. Above we have demonstrated 
that in the case of Castor the neon line ratios are influenced by 
blending with iron lines. Therefore, a more credible estimate 
for the density is derived from the spectrum after subtraction 
of the simulated iron contribution. In Table 0] we confront the 
results obtained directly from the data and after taking account 
of the simulation. It turns out that the temperature and den- 



To examine the X-ray variability of Castor A and Castor B we 
extracted their lightcurves using (A) the Chandra LETGS ze- 
roth order signal, and (B) the XMM-Newton data collected by 
Ihe EPIC pn. 

With EPIC the stars are not resolved, so that only a 
combined lightcurve of Castor AB is obtained. The EPIC pn 
lightcurve (Fig.|6l top panel) displays strong variability includ- 
ing four flares within the observing time of 50 ksec. The high 
signal-to-noise of EPIC pn allows to estimate characteristic pa- 
rameters for these flares: Rise and decay times, peak luminos- 
ity, and energy content are summarized in Table|5] All flares are 
of short-duration (less than 1 h) and of similar strength. While 
the rise takes place linearly, the decay phase follows an ex- 
ponential resulting from cooling. The decay times are typically 
about a factor 2 — 3 longer than the rise times. At peak emission 
the amplitude has increased by a factor 2 — 3. The energy bud- 
get of the flares is listed in columns 5-7 of Table|5l We give the 
total energy radiated in the flare time interval (as derived from 
an integration of the EPIC pn spectrum), the amount of energy 
above the quiescent energy during the flares, and the fractional 
increase of energy emitted during the flares (see Table caption 
for a more detailed description). 

The second to fourth flare observed with EPIC can be as- 
signed to either one of Castor A or B because of the simulta- 
neous coverage with Chandra (see Fig.|6land Table |5}- In or- 
der to determine the host of the first flare, which was observed 
by XMM-Newton alone, we make use of EPIC MOS images: 
we extract MOS 1 images for data subsets corresponding to the 
time of flares. A contour plot of each of these images is seen 
in Fig. For clarity we have underlaid the contours for the 
quiescent emission. Their two peaks are separated by ~ 3.2" 
outlining the positions of Castor A and B. A shift of the photon 
peak from one flare to the other is noticeable, indicating that the 
emission comes from different positions within the Castor AB 
system. For flare #2 to #4 the star that is closer to the peak of 
the flare contours is the one which is also listed in Table [S] 
i.e. the MOS image is consistent with the Chandra lightcurves. 
The contour of flare #1 suggests that this event took place on 
Castor A. This way we have identified the host of all flares. 
However, Fig. makes clear that this analysis puts us at the 
limit of the capability of EPIC. 

To trace the evolution of the coronal temperature during 
the flares we computed a hardness ratio derived from the 
counts in a soft band (S; 0.3 — 0.85 keV) and a hard band (H; 
0.85 — 2.0 keV) of EPIC pn. For this choice of energy bands 
the Lya line and the triplet of neon, which are indicators of hot 
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Fig. 6. X-ray lightcurves of Castor A and B on Sep 29/30, 2000. top - XMM-Newton EPIC pn where Castor A and B are unre- 
solved, binsize is 100s; middle and bottom - Chandra LETGS zeroth order signal of CastorB and Castor A; binsize is 400 s. 



Table 5. Characteristic parameters for flares observed by 
XMM-Newton EPIC on Castor A and CastorB. Rise and de- 
cay times are measured from 10 % of the amplitude to the peak 
and vice versa. Energies given in column 5-7 refer to the 0.1- 
2.4 keV band. The energy £^tot,F was computed by integrating 
the EPIC pn spectrum in the respective time interval. Ep de- 
notes the energy emitted in the flare alone, after subtraction of 
the average quiescent energy emitted over the duration of the 
flare. The subsequent column is the fraction of energy radiated 
due to the flare normalized to the quiescent energy. The flag 
in the last column identifies the host of the flare: Castor A or 
CastorB. 



of H/S is displayed in Fig. |8l For three out of the four flares 
the peak in the spectral hardness precedes that of the lightcurve 
indicating that the outburst is a result of heating. Flare #2 is re- 
markably hard, i.e. the relative amplitude of this event in H/S 
with respect to the other flares is larger than in Fig. |6l where 
count rates in the 0.3 — 5 keV broad band are shown. Note that 
this event is weaker in the LETGS lightcurve than expected 
from direct comparison with the EPIC lightcurve. It seems nat- 
ural to explain this difference by a combination of flare hard- 
ness, i.e. temperature, and instrumental sensitivity. However, 
emission at higher energies (E > 2 ke V) where the sensitivity 
of the LETGS drops sharply is not more prominent than for the 
other flares. 



# 


[min] 


'^dccay 

[min] 


^prcflarc 


lg-Etot,F 

[erg] 


IgSp 
[erg] 


Bp 

Eq 


Host 


1 


18.6 


29.8 


2.78 


32.57 


32.12 


0.63 


A+ 


2 


9.2 


15.8 


2.42 


32.27 


31.77 


0.55 


B* 


3 


11.9 


36.5 


3.11 


32.55 


32.02 


0.51 


A* 


4 


8.6 


11.6 


2.34 


32.47 


31.70 


0.58 


B* 



* Host assigned from simultaneous, resolved Chandra lightcurves, 

host assigned from photon centroid in MOS image combined with 
information for other flares from Chandra (see text in Sect.|6j. 



plasma, are located in the hard band. The temporal behavior 



A time-resolved study of the high-resolution RGS spec- 
trum would provide more information about temperature and 
density changes during the flares. However, due to poor statis- 
tics such an analysis does not seem promising. For Castor C 
which is brighter in X-rays we found enhanced emission at 
A < 16 A (either from continuum or unresolved lines), and 
marginal (1.5 cr) evidence for variable G and R ratio during 
flares JStelzer et al. 2002> . Similarly IGiidel et al. (2002)| ob- 
served variable emission line fluxes during a giant flare on 
Proxima Centauri. 
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Fig. 7. Contour plot of EPIC MOS 1 images at the position 
of Castor AB during the four flares (see Hghtcurve Fig. |6j. 
Underhned dotted contours represent the quiescent emission, 
and mark the position of Castor A and Castor B. Contours are 
normalized to the peak. The relative shift of the photon cen- 
ter indicates that the major part of the emission moved from 
Castor A to Castor B and vice versa. 
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Fig. 8. Time evolution of spectral hardness during the four large 
flai-es observed with XMM-Newton: Hard/Soft = counts [0.85- 
2.0keV] / counts [0.3-0.85 keV]. Different plotting symbols 
denote different phases of the Hghtcurve: circles - quiescence, 
asterisks - rise of flare, diamonds - decay of flare. The bins 
for the rise phase of all flares are defined from 10 — 50% and 
50 — 100% of the peak amplitude, and for the decay phase the 
other way round. For most flares the maximum hardness, i.e. 
highest temperature, is reached during the late rise phase of the 
hghtcurve. 



7. X-ray Activity Level 

X-ray flux, luminosity, and the ratio between X-ray and bolo- 
metric luminosity are frequently used to measure the strength 
of stellar activity. We computed the X-ray luminosity of 
Castor AB by integrating the time-averaged EPIC pn spectrum. 
For the ROSAT hwA we derived Igix = 28.96 erg/s, com- 
patible with earlier measurements of EXOSAT and Einstein 
JPallavicini et al. 1990b> . and ROSAT ( Schmittetal. 1994> . 

The Chandra data can be used to obtain an independent, 
estimate of the X-ray luminosity of each of the two binaries. 
Since the zeroth order does not contain any information on 
the spectral distribution of the counts we used PIMMS^ to 
derive the X-ray luminosities. We find IgLx.A = 28.9 erg/s 
and Igl/x.B = 28.8 erg/s in the ROSAT hand assuming a one- 
temperature Raymond-Smith model of kT — IkeV. While 
clearly too simplistic this rough estimate shows that the result 
obtained with both satellites are compatible with each other 

Based on the Chandra luminosities the lg(Lx/Lboi) 
ratio can be evaluated for all components of Castor AB. 
To compute the bolometric luminosity for the A-type pri- 
maries we made use of the published V magnitudes (e.g. 
Hiinsch et al. 1999 1 and the bolometric corrections given by 
iKenyon & Hartmann (1995)| The bolometric magnitudes of 
the secondaries were kindly made available by M. Giidel 
(priv.comm.). Assuming that the X-rays are emitted by the pri- 
maries leads to lg(Lx/iboi) — —6.3 and —6.0 for Castor A 
and B, respectively. This is higher than the canonical value of 
Lx/Lhoi ~ 10~^ measured for early-type (OB) stars that pro- 
duce X-rays in their winds. If instead the unresolved secon- 
daries are responsible for the X-ray emission Ig (ix/-^boi) is 
—3.7 and —3.5, as typical for late-type active stars. Thus the 
observed emission level suggests the low-mass companions as 
X-ray emitters. 

Recently, [Ness et al. (2002a)l investigated the dependence 
between coronal temperature (expressed in the Lya/r ratio) 
and activity level (represented by the X-ray surface flux) for 
a number of late-type stars. We add Castor AB to this sample. 
When converting to flux we assumed (i) that the secon- 
daries are the source of the X-rays, and (ii) that the observed 
X-rays are distributed over the surface of both late-type stars. 
Both assumptions are supported by our discussion above. Fig.|9l 
demonstrates that Castor AB are moderately active stars that fit 
well into the straight correlation between X-ray flux and tem- 
perature characterizing stellar coronae. 

8. Summary 

Castor is so far the only stellar system observed simultaneously 
with both Chandra and XMM-Newton. The contemporaneous 
observation allowed us to explore the X-ray emission from this 
multiple star with unprecedented detail using the complemen- 
tary capabilities of both satellites. 

The a Gem system, composed of the two spectroscopic bi- 
naries Castor A and Castor B, has been resolved for the first 
time in X-ray light. At a separation of 4" the only X-ray in- 
strument capable of resolving Castor A and B completely is 



http;//asc.harvard.edu/toolkjet/pimms.jsp 
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Fig. 9. X-ray surface flux versus Lya / r line ratio of oxygen as 
a measure for the coronal temperature. Fx was computed for 
the ROSAThwd (0.1 - 2.4keV). 



presently Chandra. The Chandra LETGS allowed us to sep- 
arate both the lightcurves and the high-resolution spectrum of 
the binaries, demonstrating that these stars behave like a typi- 
cal late-type coronal X-ray emitter: they exhibit strong variabil- 
ity including several flares, and their spectrum is dominated by 
emission Unes of O VIII, O VII, Fe XVII, and Ne IX. However, 
the LETGS spectrum is of low quality (due to the off-axis po- 
sition and X-ray faintness of Castor A and B). 

A quantitative analysis of the high-resolution spectrum of 
Castor AB was performed with the XMM-Newton RGS. Its res- 
olution is insufficient to resolve the two binaries from each 
other, but its higher sensitivity allows to use individual emis- 
sion Unes and their flux ratios for coronal diagnostics. We mea- 
sured the temperature sensitive G ratio and the density sensi- 
tive R ratio of the helium-like ions of oxygen and neon, and 
compared these ratios to calculations for collisional ionization 
equilibrium to learn more about the conditions in the coro- 
nal plasma of Castor AB. Simulating a neon-free RGS spec- 
trum based on the temperature structure inferred from the EPIC 
spectrum and subtracting this simulated spectrum from the data 
we have shown that the neon triplet is severely contaminated by 
a blend of iron Unes which affect the r and the i line. As a result 
the density derived from the Ne IX triplet must be regarded as 
quite uncertain. 

For oxygen we are able to derive more reliable plasma pa- 
rameters. For this element the electron temperature is ^ 2 MK 
and the electron density is (0.5...1) 10^" cm^"^. The two values 
for the density arise from differences of the radiation field of 
Castor A and B. However, the assumption of equal line ratios 
for both stars may have introduced an uncertainty of the same 
order, such that we consider the values cited above as a likely 
range for the density rather than two different values for the 
coronae of the two stars. 

Another temperature sensitive line ratio compares the 
strength of the Lya line of an H-like ion to the resonance line of 
the He-like ion of the same element. Representing two different 
ionization stages this ratio directly relates to coronal temper- 
ature. We showed that the coronal temperature of Castor AB 
as measured by Lya/r fits well into the correlation with X- 



ray surface flux observed for active stars. Note that the Lya/r 
measurements for the active stars we used for comparison (data 
from lNess et al. 2002at were all based on LETGS data, while 
our data point for Castor derives from an RGS observation. 
Previous studies of stars observed with both the RGS and the 
LETGS indicate a tendency of the RGS to measure slightly 
higher values for Lya/r than the LETGS cAudard et al. 200 T1 
IRaassen et al. 2 002 Stelzer et al. 2002 1. However, most of the 
RGS and LETGS measurements are compatible with each 
other within the statistical uncertainties. Therefore, we con- 
sider the respective calibration of both instruments sufficient 
for a direct comparison. 

The X-ray lightcurves of Castor A and Castor B show that 
both are strongly variable and frequently subject to flaring. 
With a combination of Chandra lightcurves and EPIC MOS 
images for individual flares we established the host of each of 
the six flares: Chandra directly resolves component A and B, 
while in the MOS image the photon center shifts during flares. 
Three events occurred on Castor A and three on Castor B. Two 
of the flares were observed only with Chandra (during the 
last part of the observation, when the XMM-Newton expo- 
sure had already finished). For these two flares the low S/N 
of Chandra does not allow an analysis of its characteristics. 
For the four flares that occurred within the observing time 
of XMM-Newton we made use of the excellent sensitivity of 
EPIC pn, and derived duration and energy output for each 
of them. The flare recurrence time can be estimated to one 
event every 3 — 5 h. Similarly strong activity was observed on 
Castor AB in an earlier XMM-Newton observation presented 
by [Giidel et al. (2001)] However, lacking simultaneous high- 
spatial resolution Chandra data some uncertainty remained as 
to the origin of the flares (Castor A or B). On basis of this 
eailier XMM-Newton observation Giid el et al. (20 01) supposed 
that the flare rate on Castor B is higher than on A, and flares de- 
cay faster than on A. The observations presented here suggest 
similar flare rates for both A and B. For the flares observed with 
XMM-Newton our quantitative analysis shows indeed that de- 
cay times are shorter for B than for A. But still better statistics 
are required to clearly establish a difference in the flaring prop- 
erties of Castor A and B. The flares on all visual components of 
the Castor system (including YYGem; see lStelzer et al. 2002> 
are of short-duration similar to the majority of X-ray flares seen 
on young stars, but in contrast to the common type of long- 
duration flares on more evolved stars such as RS CVn binaries. 
The X-ray luminosity and strength of all X-ray flares observed 
on Castor A and B so far is similar to that of dMe flare stars with 
amplitudes of a factor of 2. . .3 (see e.g. lPallavicini et al. 1990bt . 
This supports the common believe that the late-type compan- 
ions in each of the binaries are responsible for the X-ray emis- 
sion. 
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